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1 Introduction 



Next-to-leading order calculations of heavy-flavour production have been available 
for a long time |l|, ^ ^. These calculations exploit the fact that the mass of 
the heavy quark acts as an infrared cutoff on coUinear singularities, and thus the 
cross section has a power expansion in the strong coupling constant, evaluated at a 
renormalization scale near the heavy-quark mass m. This approach is appropriate 
when the heavy-quark mass is the only relevant mass scale of the problem, and it 
is bound to fail when the transverse momentum of the heavy quark is much larger 
than its mass. In fact, in this case, one cannot pinpoint a single characteristic scale 
in the problem, since all momenta between m and are equally involved. It turns 
out that, whether we choose the renormalization scale /i^ and the factorization scale 
/ip of order m or of order p^,, large logarithms of the ratio Pt/ttl arise to all orders in 
the perturbative expansion, and spoil its convergence. The logarithmic structure of 
the perturbative expansion for the inclusive transverse-momentum distribution can 
be classified in terms of the form a'^{aslog p^ / m)^ (which we call leading-logarithmic 
terms, or LL), plus terms of the form a^{aslogpT,/m)^ (next-to-leading logarithmic 
terms, or NLL), and so on. 

An early attempt to deal with this problem was described in ref. 0, where an 
order-Qfg calculation of the p^ spectrum was presented. There, in order to deal with 
the logarithmically enhanced terms, the following approach was adopted: the scales 
/Xr and fip were chosen of the order of p^, and the first neglected LL terms (that is to 
say, the terms of order a^{as logp^/m)'^) were estimated. Their value was used as an 
error on the resulting cross section. This procedure, of course, ends up giving large 
errors for very large p^, but a reasonable range of can be accessed in this way. 

Alternatively, the whole tower of LL and NLL corrections can be computed using 
the fragmentation-function formalism 0. This approach has however the drawback 
that it is essentially a "massless" formalism, in the sense that all contributions to 
the cross section that are suppressed by powers of m/p^ are not included. Although 
these corrections must be small when \ogpT^/m is large, it is however never clear at 
what values of Pt they can really be neglected. 

The purpose of this work is to construct a formalism in which the fixed-order and 
the fragmentation-function approaches are merged, in the following sense: 

• All terms of order a1 and are included exactly, including mass effects. 
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• All terms of order aja^ pr^/m and ala^ pT^/m are also included exactly. 

To be more specific, let us write schematically the result of the NLO calculation of 
the hadronic cross section as 

da 



A{m)at + B{m)at + 0{at) . (1.1) 

The explicit dependence upon the centre-of-mass energy -Ecm, Vt ^"^^ is not indi- 
cated, and Os = «s(/u)- The NLL resummed calculation is given by 

^ oo oo 

aAa^ log u/m)* + > hAa^ log u/m)* 



dpi 

+0{at{a, log/i/m)^) + 0{al x PST) , (1.2) 

where PST stands for terms suppressed, in the large-px limit, by powers of vnjp^^ 
irrespectively of further powers of logarithms and of a^. The coefficients and hi 
depend upon i?cm, Vt and /i. If /i ~ Pt, they do not contain large logarithms of the 
order of \o<gpr^lm. The only large logarithms are the ones explicitly indicated. Our 
approach combines the results of eqs. (|1.1|) and ( |1.2|) , giving 

da 



A{m)al + B{m)oLl + 



dp\ 

aj (tts log /i/m)' + hi (tts log /i/ x G{m,p^) 



oo oo 

,2 



i=2 1=1 

+C)(a^(aslog/i/m)0 + 0(0' x PST) , (1.3) 

where the function G{m,prj,) is quite arbitrary, except that it must approach 1 when 
m/pT^O, up to terms suppressed by powers of m/px- Observe that the sums now 
start from i = 2 and i = 1, respectively, in order to avoid double counting. Thus, this 
formalism contains all the information coming from the fixed-order NLO calculation, 
and from the NLL resummed calculation. The arbitrarity in the G factor arises from 
the fact that we do not know the structure of power-suppressed terms in the NLL 
resummed calculation. 

In the present work, we adopt the shortest path for obtaining the correct answer, 
making use of an already existing computation of the 0{al) cross section {fixed-order 
approach, or FO) and of already available computer codes to evaluate the resummed 
cross section in the massless limit {resummed approach, or RS). 
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In order to carry out our program, it is important that both the RS and the FO ap- 
proaches are expressed in the same renormahzation scheme. The commonly used FO 
approach uses a renormahzation and factorization scheme in which the heavy flavour 
is treated as heavy. Thus, if, for example, we are dealing with bottom production, 
we use as of 4 light flavours as our running coupling constant, and the appropriate 
structure functions should not include the bottom quark in the evolution. The RS 
approach, on the other hand, also includes the heavy flavour as an active, light degree 
of freedom. This problem can be easily overcome by a simple change of scheme in 
the FO calculation. Section || contains the details of this procedure. 

Once this is done, the FO calculation matches exactly the terms up to order 
in the resummed approach, in the limit where power-suppressed mass terms are neg- 
ligible. In order to subtract from the RS result the fixed-order terms already present 
in the FO, we must provide an approximation to the latter where only logarithmic 
mass terms are retained. We will call this "massless limit" FOMO. In the simplified 
notation of eqs. ( |1.1|) and (|1.2| ) we have 

A(m) = ao + PST, fi(m) = ai log/i/m + 60 + PST , (1.4) 

and the FOMO approximation is given by 

-f^ = aoaj + (ai log /i/m + 60) + C'fa^ x PST) . (1.5) 

Our final result will be given by 

FONLL = FO + (RS - FOMO) x G{m,p^) . (1.6) 

The notation FONLL stands for fixed-order plus next-to-leading logs. Formula (|1.6| ) 
is our practical implementation of eq. ( |1.3| ). 

An alternative approach to the problem has been given in ref. 0, using an ex- 
tension of a method developed for the computation of heavy-flavour effects in deep- 
inelastic scattering [|^. There, however, NLL terms are not fully included. We will 
comment on this approach in due time. 

The paper is organized as follows. We first review in sect. |^ the resummed 
approach, and in sect. |^ we describe the procedure to adopt in order to translate 
the FO result from a scheme with nf — 1 light flavours to a scheme with rif light 
flavours. In sect. ^ we give a few details concerning the calculation of the massless 
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limit of the FO calculation. In sect. ^ we check the matching between the FOMO and 
the RS calculation, by explicitly verifying that the difference RS-FOMO is of order 
Og. In sect. |]we examine the size of power-suppressed effects in order to understand 
at which value of mj-p^ the massless approach gives a sensible approximation to 
the massive calculation. The function G{m,prj,) will be chosen on the basis of the 
considerations given in this section. In sect. |^ we describe a simplified version of our 
FONLL calculation, in which only LL resummation is adopted. In sect, ^we describe 
an alternative choice for the factorization scheme for the fragmentation function. In 
sect. |]we describe our full result, for the case of bottom production at the Tevatron. 
Finally, in sect. IC we give our conclusions. 

When computing specific cross sections, we will always use the parton densities 
set CTEQ3M, which implements a correct treatment of the bottom parton density. 
Since this set uses a mass of 5 GeV for the bottom quark, we will also employ the 
same value for consistency. 



2 Review of the resummation formalism 

Essential ingredients for the resummation formalism are the perturbative frag- 
mentation functions for the parton i to go into the heavy quark h, Di{x, fi), where 
i runs over all the light partons (e.g. the light quarks and antiquarks, and the gluon) 
plus the heavy quark and antiquark. These fragmentation functions satisfy the stan- 
dard Altarelli-Parisi evolution equations, and their initial values at a scale fiQ of 
the order of the heavy-quark mass m are perturbatively calculable. In the modified 
minimal subtraction scheme (MS) they are given by 



Dhix,fio) = 6{1 -x) + 



ag(/xo)CF 
2tt 



1 + /, /^o 



log — - 21og(l - X) 



X \ ' 



(2.1) 



D,{x, /.o) = ^i^^(x^ + (1 - xf) log ^ (2.2) 

Di{x,iJ,o) = 0, ioTi^g,h, (2.3) 

h being the heavy quark and g the gluon. As usual, Cp = 4/3 and Tp = 1/2. 
The notation [/(x)]^ denotes the so-called +-distribution, whose integral against any 
smooth function g{x) is defined by the equation 

' g{x) [/(x)]+ dx = C {g{x) - g{l)) f{x) dx . (2.4) 
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The perturbative fragmentation functions (PFF), evolved up to any scale ~ Pt 
via the Altarelli-Parisi equations, can be used to evaluate heavy-quark cross sections 
in the large-transverse-momentum region by convoluting them with short- distance 
cross sections for massless partons [10|, |Tl|, subtracted in the MS scheme. The 
heavy quark is also treated as a massless active flavour, and therefore also appears 
in the parton distribution functions of the colliding hadrons and in the evolution of 
the strong coupling constantQ The differential cross section for the hadropro duct ion 
process 

H^{P^) + H2{P2)^h{P) + X , 
following the notations of |^, is given by 



(Pa dz dv dw 



dprj^dy S Ji-v+vw ^ Jvw/z 1 — v Jvw/zv 



X 



xFlj^{xi, i^p)FIj^{x2, fiF)Dk{z, fip) X (2.5) 



X 



1 f da'^(s,v)\ ^.^ , agf/in) / n 

- 1 ()(l-w)H — -Kij-^k{s,v,w;Hn,IJ'F) 

v\ dv J..^^ 27r 



having defined the hadron-level quantities 

with 5 = (Pi + P2)^ T = {Pi- Pf and U = {P2 - Pf. We also define 

pi = xiPi, P2 = X2P2, p = P/z (2.7) 

and 

t —u 

s = {j)i+P2?, t={pi-pf, u = {p2-pf, v = l + -, w = —- . (2.8) 

s s + t 

In terms of the momentum fractions Xi, X2 and z, it holds 

vw 1-V 
s = xiX2S, xi = , ^2 = —^, r- (2-9) 

ZVW Z[l — V) 

The a^{s,v) terms represent the leading-order massless-parton to massless-parton 
scattering cross sections, while the Kij^k{s,v,w; ^j^, fip) represent the NLO correc- 
tions, explicitly given in ref. |^. 

■^For this reason, this resummation formahsm is sometimes referred to (somewhat improperly) as 
the "massless scheme" . 
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The implementation of the resummation procedure has been performed in ref. 
for pp, in [1^ for 7p, and finally in |jl3| for 77 collisions. In all cases the results 
agree with the full massive ones (refs. 0, [|14| and respectively) in a region 
from two to four times the mass of the heavy quark. They have a smaller scale 
dependence than the fixed-order calculations at larger Pt, because the large logarithms 
originating from gluon emission and gluon splitting are resummed by the evolution 
of the PFF. In this region they are therefore more reliable (see ref. for a more 
complete discussion on this point). Because of the intrinsically massless nature of the 
resummation procedure, these results cannot be trusted when j9t approaches m. 



3 The change of scheme 



When performing the full FO massive calculation one can, according to ref. 
conveniently define two renormalization schemes that describe the same physics. One 
is the usual MS scheme, in which all flavours are treated on an equal footing. The 
other one, which we will call the decoupling scheme, is similar to the MS scheme 
except for its treatment of the divergences arising from heavy-flavour loops, which 
are subtracted at zero external momenta. In the decoupling scheme, for processes 
taking place at energies much below the heavy-quark mass, we can forget about the 
heavy flavour altogether. Thus, in QCD, one can define a standard MS scheme, as 
well as decoupling schemes in which the top, or the top and the bottom, or the top, 
the bottom and the charm are treated as heavy. It is common to refer to the cou- 
pling constant in the corresponding schemes as the 6-, 5-, 4- and 3-flavours coupling 
constant. 

We will now focus, for simplicity, on the case in which we have nif = Uf — 1 
massless flavours and a heavy one of mass m. The strong coupling constant in the 
decoupling scheme and in the standard MS scheme will be referred to as the riif and 
rif flavours couplings respectively. It turns out that the nif and Uf flavours couplings 
are identical at a renormalization scale equal to the mass of the heavy quark: 

a("')(/i,) = a("'')(/^K) + Oial) for /i^ = m . (3.1) 
Thus, using the renormalization group equation we find 

«^^/^K) = ai--\m) -b'^:-^ a! log (3.2) 

fit 
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= atH^)-ht^ocl\og^ (3.3) 

where 

r=i^%^. (3.4) 



Observe that in the highest-order terms in eqs. (|3.2| ) and ( |3.3| ) we do specify neither 
the nf (or nif) label nor the scale in as-, since the corresponding differences are of 



higher order. Taking the difference of eqs. (|3.3| ) and ( p.2|) and using eq. ( |3.1|) we 
arrive at 

«i"''H-"K) = c^tHl^.) - ^TAog^al + 0{al) . (3.5) 



Similarly structure functions for nu and rif massless flavours must match when fi^ = m 
||17|| . More specifically, they must satisfy the conditions 

Fj"f^(x,m2) = Fj"'f)(x,m2) for j ^ 
F^"f^(x,m2) = 

Fff)(x,m2) = 0, (3.6) 

where h stands for the heavy flavour. It should be emphasized that this is a property 
of the MS subtraction scheme, and it is no longer true in other schemes. Using the 
Altarelli-Parisi evolution equations together with the matching conditions given in 
eqs. (p.6[), one can easily find the appropriate relations between the parton densities 
with TT-if and nf active flavours for fi of the order of m. 

We begin with the Altarelli-Parisi equations for the parton densities with nf = 
n\{ + 1 flavours 



(91og/i2 2n ^ 



^ Ff\xlz^^)P^\z)^-^. (3.7) 



For /i of the order of m, neglecting terms of order a^, we get 

2 

ai"'^(m)log — 



Ft\x,y.)-Ft\x,m) = Ff\xlz,rn)P]:;\z)-. (3.8) 

zvr J ^ z 



Observe that the heavy-quark density at fi = m vanishes because of eqs. (p^). We 
can then exclude it, putting j h,h in the sum. For i = h {or i = h), eqs. ( p.6|) and 
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P) yield 



as '(m)log^ 



271 



(3.9) 



which shows that the heavy-quark density is of order ctg. An equation similar to 
eq. (p.8|) holds for rin flavours: 

,2 



ai"" \m) log /> 1 J 



(3.10) 



Taking the difference between eqs. ( pTSf ) and ( |3.10D , and using eqs. (|3^ ) and (|3.1| ) we 
obtain, for i h,h: 



Fl'''\x,ij)-Fp'>{x,m 



("If)/ 



ai""^ (m) log /.I 



27r 



^ F^-"\x/z,m) [P^f (z) - p!^'^^\z)\ ^ . (3.11) 

The only splitting function that depends explicitly upon the number of light flavours 
is the gluon splitting function. We have 

2Tp 



P^f{z)-P^--\z) 
which applied to eq. (|3.11D gives 



-5(1 -z), 



(3.12) 



Tp«^f^(m) log 



The final result is then 

p("_f) 



37r 



(3.13) 



K) 



0(as) 

F("«) + 0(a,^) forjVM/^), 



pint) _ p{ni{) 



1 - 



^— log — 

37r m 



(3.14) 



These equations are easily understood in the following way. Since at yU = m the 
heavy-quark density is zero, it must be only of order as if fi is near but not equal to 
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m. The gluon density is affected at order as by the presence of the heavy flavour. For 
fi > m, there is in fact some more room for the gluons to go into sea heavy flavours, 
and therefore the gluon density is shghtly diminished. Light flavours do not couple 
directly to the heavy one. They feel its effect only because of the diminished gluon 
density, since there are less gluons to go into light flavours. But the variation of the 
gluon density is itself of order ag, and the probabihty for a gluon to go into a light 
fermion is also of order as, so that the net effect is of order a^- 

Now we have all we need to translate the heavy-flavour cross section formulae 



from the decoupling scheme of ref. [IC] to the the standard MS scheme. First of 
all, the effect of the change of scheme on terms of order is of order a^, and can 
thus be neglected. On the other hand, the Born term, of order a^, will give rise 
to corrections of order a^. In the case of the qq annihilation Born term, the quark 
densities are unaffected at order as, and thus the only correction arises from the 
coupling constant. We have 



a, 



(0) 



-f. 



(nf) 



qq 



qq 



2Tp pt^' 
a«-^log — - 



37r 



+ 0{at). (3.15) 



In the case of the gg fusion, instead, we obtain symbolically 



gg 



gg 



2T^ , /ia' 

as — log — ^ 



37r 



(3.16) 



We thus summarize the results of the present section. In order to compute flxed- 
order heavy-quark cross sections at the O^a^) level, using the standard MS scheme 
with rif flavours for both the coupling and structure functions, the modifications one 
needs to apply to the partonic cross sections of refs. |]l], |[ are the following: 



Add a term — as^ log a^^ to the qq channel cross section. 
Add a term — ag^ log af'g to the gg channel cross section. 



Observe that, in case one uses /ip = /Ir, there are no corrections to the gluon channel. 
Furthermore, for any reasonable ranges of scales and masses, these corrections are 
not large. 
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In the following, we will always refer to the FO and FOMO calculations performed 
in this transformed scheme. We will thus always assume that and the parton 
densities Fj refer to ai"*^^ and Fj"'^^ 



4 Massless limit of the fixed-order calculation 

The massless limit of the fixed-order formulae (in the sense of eq. ( |1.5|) ) is obtained 
via algebraic methods from the results of ref. 0. As pointed out there, the limiting 
procedure is non-trivial. In fact, the partonic cross sections at order contain 
distributions, such as delta functions or principal value singularities. When taking 
the massless limit, new contributions to these distributions arise. We will not report 
here the rather boring details of the algebraic procedure we followed. We would like 
to remark, however, that its correctness can be easily checked in the following way. 
We compute a heavy-flavour differential cross section at a fixed px, rapidity, and 
centre-of-mass energy. We choose the renormalization and factorization scales equal 
to Pt- Under these conditions, the mass dependence of the result is confined to the 
partonic cross sections. In the massless limit approximations, the only remnants of 
mass dependence are in logarithms of the mass in the 0{aD terms. Thus, if we plot 
the FOMO cross section versus the logarithm of the mass, we get a straight line. On 
the other hand, if we plot the full FO cross section versus the logarithm of the mass, 
it should approach the FOMO result in the limit of small masses. This is in fact what 
we observe, as can be seen in figs. |l] and 0. 

From the figures, it is quite apparent that the massless limit, as well as its im- 
plementation for the calculation of hadronic cross sections, was carried out correctly. 
There are also some important observations to make. The FOMO cross section is 
smaller than the massive calculation. On general grounds, one would instead expect 
the massive calculation to be smaller, because of the reduction of phase space due to 
the presence of masses. Moreover, at a given p^, larger amounts of incoming parton 
momenta are needed if the mass terms are kept, which also reduces the cross section. 
These arguments are in fact correct: they should be applied to the absolute value of 
the cross sections rather than to their signed value. The value of the massive cross 
section does become smaller and smaller as we approach the threshold. It happens, 
however, that the massless limit cross section, which has constant slope, changes sign 
at some point, and then keeps growing in absolute value as we approach the threshold. 
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10'^ 10"' 1 1010'^ 10'' 1 10 

m (GeV) 



Figure 1: Comparison of the FO and FOMO differential cross sections as 
a function of the logarithm of the mass, at = 5 GeV. 

Notice that the FOMO approximation is quite inaccurate even at relatively small 
values of itl/pt,. For example, from both figs. |I] and |^ we notice that already at 
m/px ~ 1/3, the FO cross section is twice as large as the FOMO result. In ref. 
it was observed that for bottom production at the Tevatron, at ^ 10 — 15 GeV, 
the resummed cross section agrees quite well with the FO one (although no partic- 
ular reliability was expected for the resummed approach in that region). Following 
the results displayed in this section, we must conclude that such an agreement was 
accidental. 

As a concluding remark for this section, we wish to point out that, in principle, 
the FOMO result could also be obtained by numerical methods. It would suffice, for 
any given Pt, to evaluate the FO cross section for a fictitious mass much smaller than 

and then extrapolate linearly in the logarithm of the mass until the real m, as 
can be seen from the aforementioned plots. Needless to say, such a method would 
however be extremely cumbersome, slower and less accurate than the analytical one 
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VS = ISOOGeV, CTEQ3M 
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- qg channels only 



massless limit 
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Figure 2: Comparison of the FO and FOMO differential cross sections as 
a function of the logarithm of the mass, at px = 50 GeV. 



we have followed. 



5 Matching 

We now examine the matching between the resummed approach and the FOMO 
calculation. 

There are ingredients in the resummed approach that are not explicitly present in 
the FOMO calculation. These are the fragmentation functions for final-state partons 
to go into the heavy quark, and the parton density for finding a heavy quark inside 
the hadron. We take for simplicity //r = //p = \Pt\, and denote the common value 
with II. The fragmentation function for any parton to go into a heavy quark has a 
power expansion in terms of the coupling constant evaluated at the scale /i, and of 
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logarithms of /i/m: 

oo k 



Dj{x, = 5Z (^f'^\^^ ^) W — "s (/^) , (5-1) 



m 
fe=0 z=o 



that can be obtained by solving the evolution equation for the fragmentation function 
at the NLL level, with the initial conditions (p.l|) - (p.3|) . Similarly, the parton density 
for finding the heavy flavour in the hadron can be expanded in the form 

oo k 

F,{x, ^,m) = J2Yl ^'(/^)) - "'(/^) • (5-2) 

k=0 1=0 

With Fi{fi) in the argument of the coefficients, we mean that the coefficients have a 
complicated functional dependence upon the parton densities evaluated at the scale 



/i. The existence of formal expansions of the form (|5.1| ) and (|5.2| ) can be easily 
proved, by writing the Altarelli-Parisi equations in integral form, and then solving 
them iteratively. We present a more detailed argument in Appendix A. 

Once eqs. dSTll) and { f).2\) are formally substituted in the RS cross section formula, 
this formula itself becomes a power expansion of the form of eq. (|1.2|) , with the 
coefficients that depend (functionally) upon the structure functions for light partons, 
in the nf-flavours scheme, evaluated at the scale fi. The FOMO calculation has an 
expansion of the same form (truncated to order af) with coefficients that are also 
dependent upon the same light-parton structure functions evaluated at the scale 
/i. Thus, because of the next-to-leading logarithmic accuracy of the resummed cross 
section, the terms up to the order will match exactly with the FOMO calculation. 

This property of the resummed calculation is built into the resummation formalism 
of ref. In this work, we have explicitly verified this matching by comparing the 
numerical value of the FOMO and RS approaches in the limit — > 0. 

In order to do this, we need to vary the coupling constant towards very small 
values, in order to check if the difference of the two approaches is indeed of order a^. 
We thus need an explicit expression for the fragmentation functions and for the heavy- 
flavour parton density in terms of as(/i). For our purpose, it is, however, sufficient to 
have an approximation for these quantities that is valid at order Og. These are easily 

■*We observe that this property of the FOMO calculation is only valid in the modified scheme 
described in section ^ If we had used the standard scheme for the fixed-order calculation, the 
structure functions and the coupling as appearing there would be those with rif — 1 flavours. 
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obtained. For the heavy-flavour parton density we use 



Fg{x/z,fi)Phgiz 



dz 



(5.3) 



which equals the corresponding expression in eq. (|3.9|) up to terms of higher order in 
ctg. For the MS fragmentation functions we instead simply use the initial conditions 
eqs. ( |2.1| , |2.2| ) evaluated at the large scale ^ = p^. and without evolution. This is 
the first term in the expansion of eq. (5J.). We thus perform the RS computation 



with these approximate parton densities and fragmentation functions; we call this 
approximation RSA. In comparing RSA and FOMO we vary from its true value 
at the given to 1/3, 1/10 and 1/30 of it. The results are collected in fig. |[ 
The interpretation of the figure is quite obvious. The 0{a1) result in the FOMO 
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m (GeV) 



a, aJ3 





0.22 



0.21 



0.2 



0.19 



0.18 



a, a./30 



10 10 

Figure 3: Matching of the FOMO and RS computations for small Og, for 
p-r = 40 GeV and y = 0. The colliding hadrons are pp at ^/S = 1.8 TeV. 
The triangles represent the RS computation, the full circles the RSA, and 
the empty circles the FOMO. 



calculation does not depend upon the mass of the heavy quark, since it does not 
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contain any logarithm. The slope of the FOMO result is only due to 0{a^) terms, 
as in figs. |l] and 0. In fact, we see that the relative slope of the result decreases 
roughly like ctg, and the relative difference between the RSA and the FOMO result 
also decreases as a^. This demonstrates that the difference between the RSA and 
FOMO results is of order aj, as expected. 

In the upper-left plot of fig. ^ the full RS result (triangles) is also shown, evaluated 
with yUp = /ia = We see that, if the mass is not too small, the full RS result agrees 
reasonably well with the approximate RSA one, the two being of course coincident 
when m = p^, and therefore no evolution whatsoever is taking place. 

We observe that, at physical values of as, and for masses of the order of the 
transverse momentum, the RS cross section differs significantly from the FOMO one, 
in spite of the fact that their difference is formally of order a^. This is unfortunate, 
since it is precisely in this region ~ that we would like to see a smooth 
matching between the FO result and the FONLL one, and this can only happen if 
the FOMO and the RS calculations cancel to a high extent. 

6 Power effects in the RS and FOMO calculations 

When comparing and matching the FO, FOMO and RS approaches, there is a 
lot of arbitrariness in the way mass effects are treated. For example, we may decide 
to compare transverse-mass distributions instead of transverse momenta. These are 
equal for the FOMO and RS calculations, but differ in the FO case. The difference 
is shown in the two plots of fig. where the cross section is plotted as a function 
of the mass, keeping either mx or Pt fixed. If we focus upon the Born result, it is 
quite easy to understand the differences in the two plots. The massive result is closer 
to the massless one when the transverse mass is kept fixed, i.e. in the left plot. In 
fact, in this case, the differences come only from mass effects in the matrix elements, 
and these effects are not large. On the other hand, if the transverse momentum is 
kept fixed (right plot), as the mass grows the FO Born result requires more incoming 
momenta from the structure functions, and this suppresses the cross section. This 
mechanism is also at work, of course, in the NLO result. However, in this case, the 
FOMO approximation is bound to fail more radically, because of its linear dependence 
in logpT^/m: for large masses, it changes sign, which is clearly unphysical. Of course, it 
could be argued that the correct form of the logarithm in the massless approximation 
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Figure 4: FO versus FOMO at Born and full 0{a^) level, plotted as a func- 
tion of the mass and at fixed transverse mass (left figure) or fixed transverse 
momenta (right figure). 



should be something like log(p^ + rn?)/rn?, but the point is that the only use of the 
FOMO computation is to subtract it from the available RS result, which only contains 
logarithms of p-r/m. 

We will therefore proceed as follows. For a given transverse momentum Pt, the 
FO cross section is evaluated and combined, using eq. ( |L6D , with the FOMO and RS 
results evaluated instead at the corresponding ttit = \/v% + m? value. In this way 
the three calculations are performed at the same mx. Moreover, a central choice for 
the renormalization and factorization scales will be = = \/])% + m^, so that 
they coincide in the three calculations. 

Having done so, one has further freedom in the choice of the function (^(m,^^), 
the factor multiplying the RS — FOMO term in eq. (|1.6|) . G{m,prj,) should approach 
1 at large Pt- It could in fact be chosen equal to 1. However, as we showed in sect. |], 
the difference RS — FOMO, although of order , is abnormally large. We interpret 
this as a consequence of the fact that, when is near m, the massless approximation 
becomes completely meaningless. This is also visible from the left plot of fig. ^, where 
we can see that the FOMO result starts to deviate significantly from the full FO one 
when the mass becomes larger than about one fifth of m^. For a physical b quark, 
with a mass of 5 GeV, this means that we should suppress the RS — FOMO term for 

smaller than about 20-25 GeV. This seemingly large value is in fact not so difficult 
to justify. The dominant logarithmic effects come in fact from the flavour-excitation 
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and gluon-splitting graphs. In order for a gluon-sphtting phenomenon to be a truly 
coUinear process, and thus dominant over power-suppressed terms, the incoming gluon 
must carry a transverse mass that should be more than four times the mass of the 
heavy quark, in such a way that the produced heavy quark and antiquark are both 
relativistic. Most of this momentum should end up in the quark (the cases in which 
the momentum is evenly shared is suppressed by the luminosity) . A similar reasoning 
also applies to the flavour excitation case. 

We thus choose G{m,p^) = p\/{p\ + c^m^), and our final formula becomes 

2 

FONLL = FO + , J RS - FOMO] , (6.1) 

with c = 5, which suppresses the resummation correction RS — FOMO for mx < 5m. 

Before closing this section it is important to recall how all these manipulations 
involving power-suppressed mass terms only affect terms of order or higher. What- 
ever theoretical uncertainty might stem from them, we must live with it, as these terms 
have never been computed. 



7 An LL implementation of the matching procedure 

We will now focus on a first implementation of our formalism, in which only 
leading logarithms are resummed. Our resummed formula will have the following 
form 

FOLL = FO - FOMOLL + RSLL , (7.1) 

where for simplicity we have chosen G(m,pT) = 1- The meaning of the various 
terms is as follows: FOLL is our fixed-order plus leading- log result, FOMOLL is the 
leading logarithmic part of FOMO, and RSLL is the leading logarithmic part of the 
RS result. Thus FOMOLL contains terms proportional to a1 and log/i/m, but no 
terms proportional to alone. More specifically, in the notation of eq. ( |1.5| ) we have 

FOMOLL = aoal + oi log — al . (7.2) 

m 

The FOMOLL contribution precisely cancel the terms up to order in RSLL. Thus, 
the difference RSLL — FOMOLL begins with terms of order log^ /i/m. Therefore, we 
have an interpolating formula that reduces to the NLO result for small px, and sums 
all the leading logarithms in the large-px limit. Our FOLL results are illustrated in 
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figs. H and |[ From fig. |, which is obtained with a choice of scales /Xr = /ip = mx, 
we see that the FOMOLL and the RSLL curves coincide at = 0, corresponding 
to /i = mx = m. This is because their difference is made up of terms of the form 
Q!g (qs logyu/m)'^ for k >2, and thus vanishes for fi = m. They also depart relatively 
slowly from each other, since there are at least two powers of logarithms in their 
difference. As a consequence of this fact, the difference between the FO and the 
FOLL curves vanishes at small The FOLL curve undershoots instead the FO 
curve at larger p^? in agreement with what was estimated in ref. 0. 

There is a basic difference between the FOLL approach and the FONLL one. In 
the latter approach, it was mandatory to keep the same scales /ia and /Xp in all terms 
of the defining equation (|1.6|) . In eq. (|7.1|) , on the other hand, while it is important 
to maintain the same scale in the last two terms (in order for their difference to be 
of order a^), this scale could differ from the one used in the FO term. In fact, the 
difference RSLL — FOMOLL is a sum of terms proportional to (og log /i/m)*^, and 
their scale variation is thus of next-to-leading order (i.e. it contains an extra power of 
as or one less power of log fi/m). Thus, by a change of scale in the RSLL — FOMOLL 
difference with respect to the FO scale, one obtains an equivalent formula, up to the 
addition of next-to-leading logarithmic terms, which are not fully specified in this 
approximation. 

This observation leads to the conclusion that no improvement in the scale de- 
pendence can be observed if only LL resummation is performed, or better, if an 
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improvement is observed it must be accidental. The reader may now wonder where 
the improvement hes in this case of LL resummation, since undoubtedly we did resum 
a tower of dominant terms. One can easily convince oneself that the improvement 



lies in the fact that, for a choice of scales near p^, the large logarithms in eq. ( |7.1| ) are 
resummed, and thus one is justified to consider only scale variations for scales of the 
order of p^. This point may sound confusing, since for heavy- flavour cross sections 
at high transverse momenta the error is usually estimated by varying the scale in a 
neighbourhood of Pt- We remind the reader that this procedure is not obvious at all: 
since we are facing a two-scale problem, one should vary the scale between m and p^. 
The justification for using the upper range /z ~ Pt only was outlined in ref. 0. In 
that paper the scale was varied in the neighbourhood of p^, but subleading logarithms 
were estimated, and included in the error. Since they turned out to be small, they 
were no longer included in subsequent works. 

Having said that, we do not expect a scale improvement in the FOLL formula. 
The scale dependence of our FOLL result, displayed in fig. ^, is indeed not smaller 
than the scale dependence of the FO result. 

It is now time to comment on the results of ref. 0]. There, the NLO fixed-order 
result is fully included in the calculation, but the resummation is only performed 
with leading-logarithmic level accuracy. The authors do find a reduction in the scale 
dependence at moderate px, while at large p^ they find stronger scale dependence. 
From the considerations given in the present section, we conclude that, in fact, the 
scale compensation they observe can only be accidental. 
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8 Alternative choice of fragmentation scheme 

As pointed out in sect. the RS and FOMO results do not agree at the point 
m = (for the choice of scale fi = Pt), where we would like to see, instead, a 
cancellation. Furthermore, also the slopes of the two curves are quite different at this 
point. These two differences can easily be traced back to two kinds of terms that arise 
in the RS calculation when convoluting the NLO, MS subtracted massless kernels with 
the structure and fragmentation functions. The offset of RS with respect to FOMO is 
due to spurious terms, originating from the product of 0{a^) terms in the kernel 
cross sections with the non-logarithmic 0{as) term in the Dh fragmentation function 
initial condition of eq. These terms do not contain logarithms of pT/fn, and 

thus remain different from zero at = They are not correctly predicted by the 
RS calculation, which is accurate up to terms proportional to log(pT/'^)- The 
difference in slope is, instead, due precisely to terms of this kind. They are correctly 
given by the RS formalism, but not included in the FOMO result. 

While nothing can be done to solve this second problem, it is instead possible 
to try to eliminate the first one. This can be done by switching from the standard 
MS scheme to a new factorization scheme where the initial condition for the heavy- 
quark to heavy-quark fragmentation function, eq. (|2.1|) , does not contain the non- 
logarithmic 0{as) term and is therefore a delta-function for fio = m. We call this 
scheme the Delta scheme (DL). The transformation to go from the MS scheme to the 
DL scheme is described in detail in Appendix B. The effect of the scheme change can 
be appreciated in fig. 0, where the cross section is again plotted as a function of the 
logarithm of the mass. One can clearly see how the DL result agrees with the FOMO 
one at the point m = p^,, thereby providing the desired cancellation. However, the 
slope remains different. 

For comparison, we also plot the RSLL result. This does not contain the \og{pT/m) 
terms, and has therefore the same slope as the FOMO result. Its offset with respect 
to the FOMO result is instead due to the non-logarithmic term, which is present 
there. 
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Figure 7: Dependence of the cross section on the logarithm of the mass, 
at prj. = 40 GeV. Shown are the full FO result (empty squares), the FOMO 
massless limit (empty circles), the LL resummed result (downward- facing 
triangles), the RS result both in the MS (upward- facing triangles) and in 
the DL (full squares) schemes. 

9 Some numerical NLL results 

This section collects the final NLL results for the heavy quark cross section, as 
obtained from combining the fixed-order massive result with the resummed one. 

We begin by combining the various terms with the choice of scales /Ur = /ip = 
/i = rriT- For the sake of discussion, we now use G{m,pT) = 1, corresponding to the 
choice c = in eq. (|6.1|) . Figure || shows the full massive result (FO) compared with 
the combined ones, both in the MS and in the DL schemes. 

We can first of all appreciate the effect of switching from the MS to the DL 
scheme. In the latter the FONLL result coincides with the FO one at the = (and 
hence fi = rriT = m) point, as desired, while the result in the MS scheme presents an 
offset. This is directly related to the offset previously noticed in fig. ^ at the m = 
point. It is worth mentioning how the cancellation in the DL scheme still relies on our 
choice of factorization scale fXp = y/m^Tpf, which ensures /ip = m at Pt = 0. Any 
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other scale choice will destroy the matching at this point, giving results qualitatively 
similar to those obtained in the MS scheme. This will appear evident when we will 
present the bands obtained by varying the renormalization/factorization scales. 

We now describe the procedure we follow in order to obtain our final theoretical 
result, including the errors estimated by varying the factorization and renormaliza- 
tion scale. First of all, we fix discussed in sect. ^. Our prediction band is 

the envelope of the curves obtained by varying the renormalization and factorization 
scales in the range [mT/2, 2mT]. We also distinguish in our calculation the factoriza- 
tion scale for the structure functions from the one for the fragmentation functions, 
varying them independently. The resulting bands, for both the MS and the DL 
schemes, are shown in figs. ^ and |10. 



A few interesting observations can be made about these plots. First of all, the 
band of the combined result is much narrower in the large-pT region. This result had 
already been found with the resummed calculation of ref. 0, and is the expected 
outcome of the resummation of the large log(pT/''^) terms. Further, we notice that 
the combined bands match exactly the FO ones at This feature is due to 

the suppression factor G(m,pT), which annihilates completely the RS — FOMO term 
at this point. Without this suppression the band would be very large in this region, 
owing to the uncontrolled behaviour of the massless approximation at small p^. Last, 
the combined band lies in the intermediate-pT region, a little above the FO one, and 
is about equally wide. The maximum increase is of the order of 10%-20% on the 
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Figure 9: Combined cross section, both in the MS and in the DL schemes. 
Shown are the bands obtained by varying independently the renormahza- 
tion and factorization scales. A suppression factor of the small p-j- region, 
according to eq. (|6.1[), with c = 5, is included. 



NLO full massive result. We also notice that the DL and MS schemes give about 
the same result: no significant differences can be observed. 

The Tevatron results for b production are usually given in terms of a cross section 
with a cut Pt > P™"- order to give a rough estimate of the effect of resummation 
for this quantity, relative to the NLO calculation, we plot it in fig. 11. Since this 
is only included as an indication, and since the cross section is reasonably uniform 
in rapidity in the central region, we do not perform the y integration in the range 

\y\ < 0.5. 

In summary, we observe that our resummation procedure indicates the presence 
of a small enhancement in the intermediate region, followed by a reduction of the 
cross section (and of the uncertainty band) at larger p^,. The question we should 
now ask is how reliable this enhancement is. Of course, if we used smaller values 
for the parameter c, the cross section at moderate p^ would increase, thus helping 
in explaining the discrepancy between theory and the Tevatron data. However, the 
considerations given in sect. ^ do indicate that c = 5 is the conservative choice we 
should make to obtain a reliable result. 
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10 Conclusions 



In summary, we have defined a procedure for improving in a systematic way 
tlie calculation of the transverse-momentum heavy-flavour spectrum. Using already 
available calculations, this procedure has been applied to obtain the spectrum of 
heavy-flavour production with NLO accuracy in ag, and with NLL resummation of 
logarithms oi Pt^/itl. 

This calculation confirms essentially the results of refs. and 0, to which it 
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reduces in the small- and large-pT limit, respectively. In the intermediate transverse 
momentum region, when is roughly between two times and five times the heavy 
quark mass, we find a slight enhancement of the cross section with respect to previous 
results. 

We can understand this enhancement in two ways. In relation to the calculation 
of ref. 1^, we find that mass corrections at order and a^, not included there, are 
large and positive in this region. In relation to the calculation of ref. [0], instead, 
we find large resummation corrections, arising at the next-to-leading order level. For 
example, large contributions arise from the notoriously large corrections to the parton 
cross sections for the gg-^g + X process, convoluted with the gluon fragmentation 
into heavy quark. These contributions are of order (if one accounts for the fact 
that the gluon fragmentation into a heavy quark is of order as), and therefore are not 
included in the calculation of ref. 0. 

In our calculation, we had the freedom to suppress resummation effects of or- 
der higher than in the small-transverse-momentum region. We included such a 
suppression, based upon the comparison of the full massive- and massless-limit cal- 
culations at order a^. Since we found that the massless limit result is not a good 
approximation to the massive one for ^ 5m, our suppression factor begins to act 
precisely in that region. 

In the present work, we have only developed the basic formalism for a reliable com- 
putation of the heavy-quark p^, spectrum. The enhancement we found goes in the 
direction of favouring a better comparison with the Tevatron data on b production. 
Other ingredients should, however, be included in a full phenomenological analy- 
sis. A non-perturbative fragmentation-function contribution, for example, should be 
extracted from e~^e~ data (see for example [^, |T^) and carefully included in the 
computation. Similarly, theoretical uncertainties arising from the error in the deter- 
mination of the strong coupling constant, of the structure functions and of the 6-quark 
mass should be fully assessed in order to give a reliable theoretical prediction. We 
will address these problems in future work. 
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Appendix A: Perturbative expansion of the heavy-quark par- 
ton density 

We now prove the existence of expansions of the form of eqs. (|5.1| ) and ( p.2| ). We 
consider the Mellin transforms of the fragmentation functions and structure functions. 
We will use the notation 

A(/xo)= [ — a;^A(x,/xo), (A.l) 
J X 

and a similar one for the structure functions. Thus, every time that the x dependence 
is not given explicitly we refer to the Mellin transforms of the corresponding quantity. 

It is well-known that fragmentation functions at a scale yU are given in terms of 
the fragmentation functions at a scale fiQ (in matrix notation) by an equation of the 
form 

D{ii) = (l + M(a,(/i),log/i//io))/^M, (A.2) 

where 1 stands for the identity matrix, and M is a matrix whose elements are power 
expansions in as{fi) and log fi/fio, starting at order as{fi), and having in all their 
terms no more powers of log/i/yUo than powers of ctg. We choose /io = m, replace 
D{m) with its power expansion in terms of as(m), which contains no logarithms, 
and replace as{fn) with its power expansion in terms of as(/i), which does contain 
logarithms, but always with a smaller power than the corresponding power of a;s(/i). 
With these operations we put -D(/i) in the form of eq. 

For eq. (|5.2| ), we begin with a similar equation for structure functions 

F(/io) = (l + M(as(yu),log/i//io))i^(/i) , (A.3) 
choose = and separate the h and h and the light components / 

Fh/-h{m) = Ff,/-M + Mf,/hj{as{fi),\ogfi/m)Fj{^) , 

Fi{m) = Fi{ij) + Mij{a,{fi),\ogfi/m)F,{fi) . (A.4) 

Observe that the index j runs also on h and h. The left-hand side of the first two 
equations can be given as a power expansion in as{m), with coefficients proportional 
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to the Fi{m). They in fact start at order in the MS scheme. The Fi{m) can 
be obtained from the remaining equations. Furthermore, asim) can be expanded in 
powers of a^ifi) as before. The first two equations can then be solved for Ff^/j^{fi) 
(since the coefficients in front of the F^jj^i^^) terms begin with 1, and can therefore 
be inverted). The expressions for the Ff^ij^{^) components will be linear functions of 
the Fi{^), with coefficients that are power expansions in as{fi) and in log/i/m, with 
no more powers of logarithms than powers of as{fi), as in eq. { f)-'2\ ). 



Appendix B: The DL scheme 



We now discuss the change of scheme from the MS to the DL scheme. This will be 
given, for the fragmentation functions and for the next-to-leading evolution kernels, 
in terms of Mellin transforms, defined according to eq. ( |A.! 



Using the notation of [^, we can describe the evolution in terms of the variable 

(B.l) 



instead of /x. Calling i?(/i,/io) the evolution matrix from /io to /x (in space), we 
have 



dE{^, /ip) 
dt 



R = P- ^Pa . (B.2) 
The fragmentation functions in the MS scheme are given at the low scale by eqs. 



( |2.1| )- (|2^ ), which we rewrite here as 

a(/io) 



where 



di = 



1 + 



27r 



P,,log4 + ^d, + Oia 



27r 



27r 

2^ 



0{a 



(B.3) 



dx 

X 



— X Cp 



l+x^ 

1 — X 



(21og(l-x) + l) 
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gq 



99 



dx 



X 



1 + x^ 

1 — X 



We want a scheme D where 



I ^x''T,[{l-xY + x'] . 



DM = 1 + 



(B.4) 



(B.5) 



and the other components of the fragmentation function remain identical. 
The transformations 
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have precisely this property, that is to say they translate eqs. (p.3|) into 



(B.6) 



27r 



(B.7) 



Having changed the initial conditions, we now need to find the modifications to 
the evolution equations and to the short- distance cross sections. Defining the matrix 



-5ijdi hr g, C, 



ig ^gi 



0, c,, = o. 



we have 



D 



In 



(B.8) 
(B.9) 



where we imply vector and matrix notation when we omit the indeces. 
Since physical quantities should remain unchanged, we must have 

To satisfy the above equality, we define 



(B.IO) 



1 + 



a I 



2n 



2tx 



Cj i?(/i,/io) 
C 



2tx 



(B.ll) 
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We now derive the evolution equation for E. We use the obvious observation that 
(1 — a{^) /2t[ C)E obeys the same evolution equation as E 



di 
Since 



we get 



27r 



C E{n,no) 



27C 



27r 



C E{n,no) 



dt 



bo , . 



(B.12) 
(B.13) 



bp ajfi) 

2 27T 



CE{fx,fxo) + 1 



2n 
2tx 



C 
R 



dt 

l-^C]E 
2-K 



(B.14) 



Multiplying both sides on the left by (1 + a(/i)/27rC), and neglecting higher-order 
terms consistently, we get 



d-E(/i,/io) 

dt 



2n 



R + CPo- PnC - 



£■(/!, /io) . 



[B.15) 



As a last point, we need to specialize our result to the standard cases in which one 
separates singlet and non-singlet component in the evolution equation. It is clear 
that for the non-singlet components we have Cns = —di, and the commutator term 
is zero. For the gluon and singlet components, the Cs matrix has the form 



^99 


^gq 










^q9 


Cqq 







—di 



(B.16) 



Next we describe how to implement the scheme change in the short- distance 
partonic kernel cross sections, written in the notation of 0. The Born cross section 
has the form (see also eq. ( p.5|) ) 

■'1-(1-V)/X3 



dx-^ 



where 



VW/x-i 

VW 



1 - V 



Xi 



X2 



VXs [1 

Observe that the limits on v imply a limit on x^ 



V)X3 



VW 1 - V 

< 1 >X3>l-V + VW 



X3 



X3 



(B.17) 



[B.18) 



(B.19) 
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We now introduce a scheme transformation for the fragmentation function 

D{xs) = Dixs) + Jdydz 6ix3 - yz) p{y)b{z) . (B.20) 

The second term of eq. ( |B.2(J| ), combined with the Born term in eq. ( [B.17| ), gives rise 
to a correction to the next-to-leading contribution. We wish to cast this correction 
in the same form as the next-to-leading term in eq. ( p.5|) , so as to implement it easily 
in the numerical code. 

The following change of variables 

^' = ^7T v' = l-y{l-v), (B.21) 

vy + [l- y) 

with the inverse 



v'w' 



v'w' + I — v' 

maps the unit square onto the unit square, thus the integration range < y < 1, < 
f < 1 is equivalent to the range 0<w'<l,0<f'<l. The determinant is 

v' 



and 



The limits on v imply 



dv dy = dv' dw' — (B.23) 



1-V 

xi = , X2 = — . B.24 

V W Z [1 — v'jz 



1-V + VW < z < 1 , 

VW , . 1-^ 
< v' < 1 



z 
VW 



< w' < 1 . (B.25) 



In terms of the new variables the cross section becomes 

-1 rl-{l-V)/xa 



ab = dxs / dvFQ{v, X3)C{xi, X2)D{x3)/xl + 5a, (B.26) 

Ji-v+vw Jvw/x3 



where 



5a = I dz dv' / dw' 

n-v+vw Jvw/z JvW/(zv') 



Fo{v, yz) C{xi, X2) p{y) D{z)/z^ . (B.27) 



{v'w' + 1 - v')^ 
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The structure of this equation coincides with that of the next-to-leading term in 
eq. (|2.5| ), as desired. 

Since p{y) is a distribution at y = 1, it is convenient to rewrite the above formula 

as 

fi 

dw' p{y) 



5a = j dz dv'< / 

'l-V+VW JvW/z { JvW/{zv') 



V' 



(v'w 



— — ^^o(^^, yz) C{xi, X2) - v'Fo{v', z) X2) 

' + 1 — V y 

l'l~(l~V)/z 

+ I dz dv'v'Fo{v',z)C{x[,X2)D{z)/z'^ 

'i-v+vw Jvw/z 

X / dw' p{y) 

JvW/{zv') 



b{z)/z^ 



where 



In the DL scheme 



p{y) 



aCf_ 
27r 



^1/2 ~ ^l/2\w'=l 



-l-21og(l-i/)) 



Using the fact that 



/ p{y) 

Jo 



dy = 



we have 
»i 



dw' p{y) 



'VW/(zv') 

where 



VW/(zv') 



piv'w' + 1 — f ') dw' 



j-l-v'+VW/z 



V' 



(B.28) 
(B.29) 

(B.30) 
(B.31) 

p{y) dy , 

(B.32) 



j\{y) dy=^ (log(l - y){y' + 2y-l)-2y + 21og2(l - y)) . (B.33) 
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